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ABSTRACT 



The design and construction of a data recording and processing 
system is described which permits simultaneous measurements from three 
proton free-precession magnetometers to be written onto tape at two 
second intervals. Detailed circuit diagrams and a computer program 
used for de-bugging the magnetic tape output are presented. 
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I . INTRODUCTION 



Since free nuclear precession in a weak magnetic field was first 
observed by M. Packard and R. Varian in 1953 [1] various magnetometers 
have been developed and constructed to measure magnitude, direction, or 
variation of the geomagnetic field. 

A widely used magnetometer is the Proton Free-Precess ion Magnetometer. 
Its advantages are low cost, simplicity, ruggedness, and easy portability. 
In addition it needs no calibration and is virtually unaffected by changes 
in orientation and temperature. However, its sensitivity is much lower 
than that of the more sophisticated rubidium, helium and cesium devices. 

The application of magnetometers in geomagnetic gradiometers is 
based on the assumption that the earth’s magnetic field is invariant 
over short time periods (see i.e. Refs. 2-8). Witham and Niblett [9] 
have analyzed the errors introduced into aeromagnetic survey by geo- 
magnetic time variations. They found that the rms differences in time 
variations are proportional to the rms level of magnetic activity and 
to the separation of two stations 50 feet and 26 and 94 miles apart. 

Duffus et. al. [10] have determined the short-range spatial coherence 
of geomagnetic micropulsations for stations 10, 50, 300, 950, and 1800 
kilometers apart. The authors concluded that coherence was high only 
for moderate signal levels. Similar work was done by Jacobs and Sinno 
[11] and Wright [12]. So far, however, it is not known that detailed 
analysis has been carried out to show to what extent the assumption of 
short-time invariance of the geomagnetic field is valid for very short 
station separations. For this reason a solid state data recording and 
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processing system with three identical, individually controlled and 
monitored proton free -precess ion magnetometers was designed and con- 
structed for the study of geomagnetic field coherence. The system allows 
simultaneous measurements at two second intervals from three magnetometers 
to be written onto magnetic tape for processing on digital computer and 
data analysis. 
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II. BRIEF INTRODUCTION TO THE PRINCIPLE OF OPERATION 
OF THE PROTON FREE -PRECESSION MAGNETOMETER 



The Proton Free-Precess ion Magnetometer basically consists of a 
polarizing and sensing coil surrounding a proton sample such as kerosene. 
A dc polarizing current of two to three amperes applied to the coil with 
its axis oriented approximately normal to the geomagnetic field lines 
will produce a magnetic field in the sample which is about 200 times 
stronger than that of the earth. This field will cause the protons, 
which can be considered as tiny bar magnets, to align their longitudinal 
spin axes parallel to the magnetic field lines of the coil. When the 
polarizing current is turned off, the protons tend to realign themselves 
with their spin axes parallel to the remaining geomagnetic field lines. 

During the process of realignment the protons precess coherently 
about the earth's magnetic field vector, H^, at a frequency directly 
proportional to the magnitude, |h^|, of the geomagnetic field. This 
frequency of the proton precession can be used to determine the absolute 
geomagnetic field strength, if no other magnetic field than that of the 
earth acts on the protons, by: 

“l = 2tt£ l = 7 p H i (1) 

where : 

f^ = Larmor, or precession frequency of the proton, 

7 ^ = gyromagnetic ratio of the proton, 

= proton magnetic moment/proton angular momentum. 

The value of 7 used in this paper is: 

P 

7 ^ = (.267521 + .000002) gamma ^sec ^ 
as found by Thomas [13] for mineral oil. 
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The precession signal as induced in the sensing coil is an exponenti- 
ally decaying sinusoidal voltage in the microvolt range and is of the 
form : 

v(t) = Ae a<= sin + cp(t)^ 

where : 

A = amplitude factor 
a = attenuation constant 
cp = phase? shift. 

For an earth f s magnetic field strength of about 50,000 gammas 
(1 gamma = 10 ^ Gauss), f is of the order of 2 kHz, since for protons 

i_i 

in kerosene (mineral oil) 

7 _ 

f_ = 4 . 3 kHz/Gauss. 

L zrr 

The technique used to measure the geomagnetic field strength is a cycle- 
counting method rather than a frequency determining method as will be 
shown below. 
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III. DATA RECORDING AND PROCESSING SYSTEM 



A. GENERAL 

Based on earlier work of Anderson [14] and Hansen [15] a data 
processing and recording system was designed and constructed which 
allows simultaneous readings from three proton free-precess ion 
magnetometers (Varian sensing head X 49 - 406A) to be processed and 
written onto magnetic tape for convenient data analysis. 

Although the old system was still semi-operable, considerable time 
and effort were expended to significantly update the equipment, resulting 
in improved reliability and much higher resolution. All surplus trans- 
istors, which were demonstrating high failure rates, and all homemade 
circuit boards were replaced. This made it necessary not only to re- 
design most of the circuitry, but also to replace many system components. 

A major interfacing problem developed when attempts were made to 
systematize old and new components. This resulted in a progressive 
replacement and renewal process. In the end, the only equipment re- 
tained from the old system were the magnetometer heads and the incre- 
mental tape recorder. 

The exclusive use of Flip Chip Modules (Digital Equipment Corporation) 
in the logic circuitry provided the answer to the problem of reliable 
data processing and recording in a noisy environment. 

Power requirements for the system outlined below are +10 volts at 
.5 ampere, -15 volts at 3 amperes, and -20 volts at 10 amperes. The 
size of the complete unit (except magnetometer heads, amplifiers, and 
incremental tape recorder) is 40 x 40 x 20 cm. 
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B. SYSTEM DESCRIPTION 



Figure 1 shows a simplified block diagram of the data recording and 
processing system. 

100 nsec pulses with a repetition rate of one every two seconds 
trigger three -20 volt/10 ampere drivers. The drivers simultaneously 
energize the polarizing coils of the magnetometers for a predetermined 
time (.4 to 1.3 second). Mercury wetted reed relays are used to prevent 
overloading the pre -amplifiers during the polarizing period. 

i 

After the polarizing field has collapsed and all transients from 
the magnetometer coil have died out the amplified proton precession 
signal with a usable duration of about one second is applied to the 
counting circuit the output of which is then written onto magnetic tape. 

1 . Timing Circuit 

i ■* 

A block diagram of the timing circuit is shown in Figure 2. 

Every two seconds the .5 Hz timer resets all counters and 
starts two monostable multivibrators (one-shots). One of the multi- 
vibrators produces a 800 msec pulse used to activate three relay 
drivers which open the reed relays. The 10-30 msec output pulse from 
the second multivibrator initiates another, 750 msec, one-shot which 
controls the polarizing current to the magnetometer coils. The relays 
are thus opened 10-30 msec before the polarizing period starts and 
connect the sensing coils to the pre -amplifiers only 20-40 msec after 
the polarizing current is turned off. An additional variable mono- 
stable multivibrator permits delay of the start of the counting cycle 
further if the transients from the collapsing field in the magnetometer 
coils extend beyond 20-40 msec. 
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Figure 2: Timing Circuit, Block Diagram 



Timing of all multivibrators is variable externally such that 
the precession signal amplitude and usable duration may be maximized 
according to experimental requirements. 

2 . Counting Circuit 

As was pointed out above the geomagnetic field strength is 
determined by a cycle-counting method. 

The usable duration of the proton precession signal of only 
about one second limits the time available to measure its frequency 
very accurately. For a geomagnetic field strength of about 50,000 
gammas (the actual value in Monterey is about 50,950 gammas as found 
in Ref. 15) the signal frequency from (1) is approximately 2150 Hz, 

To detect a field variation of one gamma will thus require that the 
precession frequency be measured with an accuracy of + . 04 Hz , a value 
which can hardly be achieved over a measurement time of around one 
second. This means that the error in field strength measurements due 
to the least cycle count would be + 25 gammas, which cannot be tolerated. 
The most commonly used technique to improve the measurement accuracy is 
to count the number of high frequency pulses from a very stable crystal 
oscillator clock that occur during the time required to detect a fixed 
number of cycles (1024 in the system described) of the precession 
signal. For 1024 cycles of the precession signal with a period of 
1/2150, or .465 msec and a 300 kHz clock frequency, the accuracy of 
geomagnetic field strength measurements will be approximately one part 
in 1.4 x 10^, in other words, geomagnetic field variations of + .35 
gamma or less can be detected. 

A block diagram of the counting circuit is presented in Figure 
3. 
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Figure 3: Counting Circuit, Block Diagram 



20-40 msec after the polarizing current to the magnetometer 
coils has been turned off the amplified proton precession signal is 
applied to a Schmitt trigger which produces constant amplitude square 
pulses at the same frequency. Switching threshold levels of the Schmitt 
trigger were set to -2.2 volts and -0.8 volt to prevent, as much as 
possible, triggering by noise. Any input signal between + 10 volts will 
not cause damage to the circuitry. The output of the Schmitt trigger 
is fed to an AND-gate which is enabled by control flip-flop 1 after all 
transients from the collapsing polarizing field have died out. The 
first pulse into the ten-bit slow counter sets control flip-flop 2 
which in turn enables a second AND-gate. This AND-gate allows the 
18-bit fast counter to count 100 nsec pulses from a 300 kHz crystal 
clock with an accuracy exceeding .01$. After 2^ cycles of the 
precession signal the slow counter resets control flip-flops 1 and 2, 
and the fast counter will hold the number of 300 kHz pulses counted 
for processing and writing onto tape. 

If control flip-flop 2 is in its high state for 1024 cycles of 
the return signal from the magnetometer head the count of the fast 
counter, M, can be found from: 




where : 



f c = clock frequency = 300 kHz. 



Combining (1) and (2) ; 



o O 10 -P 

2tt 2 f 



„ M , — _ 



gammas 



,8 



,, , N - 72 x 10 
H(7) = gammas 



( 3 ) 
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Neglecting the propagation delay in the slow and fast counters which is 
.7 and 1.3 pisec , respectively, and the gating speed of only a few 
hundred nsec it is apparent from (3) that the accuracy of geomagnetic 
field strength measurement's depends largely upon M, hence the + last 
count error of the fast counter and the accuracy of the crystal clock. 

The ultimate limit on the accuracy, as derived in Ref. 15 is set by the 
signal- to-noise ratio of the precession signal which has to be at least 
10 : 1 at the end of the last count in the slow counter. 

' 3 • Interface and Gating System 

Figure 4 is a block diagram of the interface and gating system. 

At the end of every counting cycle the slow counter of Magneto- 
meter 1 starts a "step-write" cycle and 

- sets the first bit of a ten-bit shift register to its high state and 
' bits two through ten to their low states, 

- fires two monostable multivibrators, 

- starts a record gap timer. 

The first nine bits of the shift register are used to sequentially enable 
six parallel AND-gates each so that the contents of the three fast 
counters are written onto tape in the order: bits 13-18 Magnetometer 1, 

bits 7-12 Magnetometer 1, bits 1-6 Magnetometer 1, bits 13-18 Magneto- 
meter 2, etc. The tenth bit of the shift register terminates the "step- 
write" cycle. 

The 10 msec one-shot steps the magnetic tape .005 inch on each 
trigger pulse received assuring that the tape is always in position 
when the next data character is up in the gating system. 

The 3 msec one-shot, together with a second 10 msec one-shot, 
generates the shift pulses to the shift register. In this manner the 
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Figure U: Interface and Gating System, Block Diagram 



shift pulses are delayed with respect to the ,, step u pulses by 3 msec 
which allows data to be written onto tape only after the tape is 
pos itioned . 

The record gap timer is a ten-bit up-counter which advances the 
tape .75 inch once every 2^ "s tep-write" cycles. 

Included in Appendix A are detailed circuit diagrams for the 
complete logic system described above. 

For further information on the logic modules used in the con- 
struction the interested reader is referred to the 1968 Logic Handbook 
issued by Digital Equipment Corporation, Maynard, Massachusetts. 

C. DATA OUTPUT 

An incremental tape recorder [15] was employed to write data on a 
standard seven-track, 1/2 inch magnetic tape with a density of 200 bits 
per inch. The tapes contained the measurement information in the fol- 
lowing format: three six-bit characters for Magnetometer 1, three 

characters for Magnetometer 3, three characters for Magnetometer 3, 

followed by a 11 (all ones) character for identification, or, a total 
o 

of ten characters per data set once every two seconds, with record gaps 
every 2^ data sets. A daily record of geomagnetic field strength 
measurements thus consists of 43,200 readings per magnetometer, or 
432 ,000 characters . 

It was found that due to noise on the 60 Hz power line the incre- 
mental tape unit was triggered from time to time accidentally causing 
more than one 77 Q characters to be written at the end of a data set. 

o 

This, however, did not affect the quality of the records and could be 
taken care of by proper programming. 
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Because of the vast amount of data per daily record it was decided 
to use the SDS 9300 computer for intermediate processing and debugging 
only, and that future coherence analysis should be done on the IBM/360 
computer ♦ 

Processing and debugging of the data on the SDS 9300 used a priori 
knowledge about the magnetometer readings . The average earth magnetic 
field strength at Monterey is about 50,950 gammas [14] with daily vari- 
ations claimed not to exceed + 100 gammas. From (3) this field strength 
is related to the number of 300 kHz pulses counted in the fast counters 
during a measurement cycle, M, by: 

2TT x 2 10 f c 1 

M = ~ p X H( 7 )’ 



so that M should be restricted to a range of approximately: 



424000 o < M < 425100 . 

O — — o 

Each data set will therefore start with an 42 character followed by 

o 

another eight information characters and an indetermined number of 77 
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characters. An SDS 9300 word holds four bytes and data will be stored 
in buffer in the form: 



| < — MAG ] 


L — ^1 


l < 


-MAG 2— 


>| 


|« 


-MAG 


3— 


> 1 






42 8 


X 

X 

00 


XX 8 


00 

CM 




^8 


^8 


42 8 


^8 




^8 


oo 

r*- 


?? 


?? 



indetermined number 77 

o 

Ideally (each data set being ten bytes only) each five computer words 
will consist of two readings from each magnetometer including the 
identifying character . 

A computer subroutine "MAGDATA" (see Appendix B) was written using 
the following philosophy: 
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- search for the first 42 character 

o 

- take this and the next two characters and store them as 0042XXXX , or 

o 

v the first data word 

- make the next three characters data word 2 and store 

- make the following three characters data word 3 and store 

- dump all 77 characters following the third data word- 

o 

- look' for the next 42 character following the last '77 character 

o o 

dumped , • 

- repeat the above up to record gap 

- when end of record sensed write debugged data onto new tape 

- continue until end of file. 

Since the record gap timer is simply a ten-bit up-counter, each record 

contains 1024 readings/magnetometer, or a total of 3072 data words. 

Writing record gaps every 10240 characters required 6000 o storage 

o 

locations each in the input and output buffers of the SDS 9300. An 

overall of 43 records/day were processed and debugged by this procedure 

and used not more than about one minute of computer time. 

The only disadvantage of the above method is that in case the fast 

counter of Magnetometer 1 fails to function properly and the first 

character of a data set is not 42 OJ all data obtained is useless. It 

o 

is therefore recommended that at least bits 13 through 18 of the fast 
counter of Magnetometer 1 be checked before and after each daily record 
and that a printout of the first and last records of each debugged tape 
be obtained before further processing. 
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IV. RESULTS AND CONCLUSIONS 



The data recording and processing system described in this paper 
represents, by far, not the ultimate possible in compactness, overall 
accuracy, and cost. However, a considerab le. portion of the Flip-Chip 
Modules was already available and had to be used in order to avoid 
excessive expenditures. The construction faced severe time and supply 
problems, to date i.e., the crystal clock is not yet available so that 
final conclusions as to the absolute accuracy of the system cannot be 
drawn . 

All components used proved to be very reliable, in more than 500 
hours of operation no component failures have been observed.. Contlnuoais 
operation over a period of more than 120 hours with an ideal input signal 
substituting for the precession signal and a 300 kHz clock with an 
accuracy better than +10 $ showed identical readings from all three 

fast counters with a + 1 count deviation as low as one per 100 measure- 
ments. The system was capable of following frequency changes of the 
ideal input signal down to + .02 Hz instantaneously, provided the Input 
signal exceeded the threshold levels of -2.2 and -0.8 volts and supplied 
at least 2 milliamperes current per Schmitt trigger. 

So far., actual geomagnetic field measurements could not be carried 
out because the pre -amplifiers and amplifiers built in 1964 [14] proved 
to be highly unreliable. The amplifiers showed the tendency to break 
into self-oscillations upon reception of any spurious signal from the 
magnetometer heads, such as transient noise, or relay noise. The actual 
precession signal could be observed by oscilloscope but was completely 
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masked by noise. Future incorporation of high performance operational 
amplifiers, such as the only recently marketed Fairchild jjA. 741, should 
solve the problem of signal acquisition reliably. It is suggested that 
the pre -amplifiers be placed close to the magnetometer heads in a 
magnetically clean box and that passive filtering be used immediately 
following the pre -amplifiers to assure the required s igna 1-to-noise 
ratio [15] at the input to the Schmitt triggers. A further improvement 
is possible by using separate charging and signal lines which was not 
attempted in the system outlined above. Since the transverse relaxation 
time of the protons is a direct function of the polarizing current, con- 
cern should also be given to the problem of increasing both the magnitude 

i 

and the duration of the polarizing current which is possible within the 
limits of -12 to -25 volts and 400 to 1300 msec, respectively. 

It is felt at this point that for preliminary geomagnetic field 
coherence studies the present system accuracy suffices. This accuracy 
will be limited by the + .01$ accuracy of the crystal clock and corres- 
ponds to a measurement error of + 5 gammas. For future experiments, 
particularly where the target detection problem is the major concern, 
the above error is intolerable and the implementation of a clock with 
an accuracy of at least one part in 10^ becomes unavoidable. 

To further increase the system sensitivity the clock frequency could 
be increased to 1 MHz which would necessitate extension of the fast 

counters beyond 18 bits and require that the clock accuracy be one part 

6 

in 10 . 

In conclusion the primary features of the system described here are: 

1. all components are plug-in units for easy replacement 

2. the system is rather compact 
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3* extension of the system in order to achieve greater accuracy and/or 
higher sensitivity is relatively easy 

4. the timing circuitry is very flexible, with the exception that the 
data rate cannot be increased 

5. reliability seems to be high 

6* power requirements are moderate 

7 \ data output is magnetic tape for easy computer processing and 
correlation analysis. 
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DETAILED CIRCUIT DIAGRAMS 



Listed below are the DEC logic and other modules used in the con- 
struction of the system described above, together with their basic 
characteristics. Standard voltage levels are -3 volts and ground, logic 
levels are: logical 1 = -3 volts, logical 0 = 0 volt. All one-letter 

* f j k . 

designators in the detailed circuit diagrams correspond to pin connections 
on the modules. 



MODULE 


FUNCTION 


CHARACTERISTICS 


R 107 


INVERTER 


Output standard levels -3 volts and 
ground 


R 113 


NAND/NOR GATE 


NAND gate for negative inputs, NOR 
gate for ground levels 


R 121 


NAND/NOR GATE 


2, 3, or 4 input gate, otherwise 
same as R 113 


R 131 


EXCLUSIVE OR 


Output -3 volts for both inputs 
the same, output ground for differ- 
ent input levels 


R 200 


RS FLIP-FLOP 


SET input makes 1 output go to -3 
volts, RESET input makes 1 output 
go to ground. Inputs are 100 nsec 
pulses -3 volts to ground 


R 202 


DUAL FLIP-FLOP 


Applications : up counter and shift 

register. SET, CLEAR, and SHIFT 
inputs 400 nsec standard pulses -3 
volts to ground. 0 outputs used in 
18-bit fast counters to match levels 
required for tape deck input 


R 302 


DELAY (ONE-SHOT) 


Normal level is ground level. Out- 
put goes to -3 volts for predeter- 
mined time when input is triggered 
by positive going level -3 volts to 
ground 


R 401 


VARIABLE CLOCK (TIMER) 


Output standard 100 nsec pulses -3 
volts to ground at any frequency 
between 30 Hz to 2 MHz 
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MODULE 


FUNCTION 


CHARACTERISTICS 


R 405 


CRYSTAL CLOCK 


Output 100 nsec pulses -3 volts to 
ground at specified frequency (5 
kHz to 2 MHz). Frequency constant 
within + .01$ for temperatures 
between 0° and +55° C 


R 601 


PULSE AMPLIFIER 


Output 100 or 400 nsec pulses -3 
volts to ground with standardized 
amplitude and width. Enabled by 
positive going levels, -3 volts to 
ground 


R 603 


PULSE AMPLIFIER 


Output 100 nsec pulses, otherwise 
same as R 601 


W 042 


10 AMP DRIVER 


10 amp of dc drive shared between 
4 outputs. Negative input brings 
corresponding output to external 
supply voltage level between -12 
and -23 volts 


W 051 


INDICATOR/RELAY DRIVER 


-3 volt input brings output to 
ground. External load connected 
to any voltage between 0 and -12 
volts 


W 501 


NEGATIVE INPUT CON- 
VERTER AND SCHMITT 
TRIGGER 


Ground level input produces -3 
volts output, and vice versa. 
Nominal switching thresholds be twee 
0 and -2.5 volts, minimum separatio 
.5 volt 


W 600 


NEGATIVE OUTPUT CON- 
VERTER 


Converts input levels of -3 volts 
and ground to outputs of ground and 
an externally supplied voltage 
(between -1 and -15 volts), 
respectively 
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TIMING AND COUNTING CIRCUIT 
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COMPUTER SUBROUTINE 



* SDS 9300 SUBROUTINE "MAGDATA" FOR READING RAW DATA FROM RECORD TAPE 1, 

* DEBUGGING THEM, AND REWRITING THEM SEQUENTIALLY IN BINARY ONTO TAPE 2 



A 


EQU 


5 


B 


EQU 


h 


READ 


TRT 


0,1 




CAT 


0 




BRU 


$-2 




LDA 


»0U0105777 




STA 


IN 




XMX 


IN,1 




C0PY 


(0,5) 


AGAIN 


STA 


IN,1 




STA 


QftjT,l 




BRX 


AGAIN,1 




LDB 


-IN 




LLSD 


10 




LRSB 


10 




LLSA 


5 




MRG 


E0M 




ADD 


*3 




STA 


ALCIN 




STB 


INCW 




RTB 


*o,i,U 


ALCIN 


E0M 


oiUooo 




P0T 


INCW 




CAT 


0 




BRU 


$-1 




ASC 


0 




PIN 


LASTW 




LDA 


•IN 




STA 


ADIN 




LDA 


-$UT 




STA 


AD0UT 




ERU 


SEARCH 


I NCW 


PZE 




SEARCH 


LDA 


*ADIN 




LDB 


-077000000 




SKM 


-0U2000000 




BRU 


«+2 




BRU 


SRCH1 




LDB 


*077770000 




SKM 


*0771*20000 




BRU 


1+2 




BRU 


SRCH2 




LRSB 


6 




SKM 


*000771*200 




BRU 


$+2 




BRU 


SRCH3 




LRSB 


6 



READ TAPE 1 



ADIN IS INPUT BUFFER 
ADOUT IS OUTPUT BUFFER 

SEARCH FOR FIRST DATA WORD 
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WRITE 



WRIT1 

ALC0UT 



E0M 

0UTCW 

SRCH1 



SKM =00000771*2 

BRU $+2 

BRU SRCHl* 

MP0 ADIN 

LDA LASTW 

ETR *077777 

SKU ADIN 

BRU WRITE 

SKG ADIN 

BRU WRITE 

BRU SEARCH 

TRT 0,2 WRITE TAPE 2 

CAT 0 

BRU $-2 

C0FY (0,5) 

LDP =0UT 

LLSD 10 

LRSB 10 

LLSA 5 

MRG E0M 

STA ALC0UT 

STB 0UTCW 

LDA AD0UT 

SUB =0UT 

C0PT (0,B) 

LRSD 10 

C0PY (A,B),(B,A) 

ADM 0UTCW 

C0PY (B,A) 

ADM ALC0UT 

LDX =00177770,2 

WTB #0,2,1* 

PZE 

P0T 0UTCW 

CAT 0 

BRU $-1 

CET 0 ERROR TEST TAPE 2 

BRU Ul 

BRU READ 

SRB #0,2,1* 

EXU ALC0UT 

P0T 0UTCW 

CAT 0 

BRU $-1 

BRX WRIT1,2 

EFT #0,2,1* 

EXU ALC0UT 

P0T 0UTCW 

CAT 0 

BRU $-1 

BRU WRIT1-1 

EjAl 011*000 

PZE 

C0PY (0,1*) 
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SRCH2 



SRCH3 



LRSD 


6 


STA 


•WAD0UT 


MP0 


AD0UT 


MP0 


ADIN 


C0PY 


(0,5) 


LLSD 


6 


LDB 


*ADIN 


LLSD 


12 


STA 


*AD0UT 


MP0 


AD0UT 


C0PY 


(B,A) 


MP0 


ADIN 


LDB 


■KADIN 


LRSA 


12 


LLSD 


6 


STA 


*AD0UT 


MP0 


AD0UT 


LDA 


LASTW 


ETR 


-077777 


SKU 


ADIN 


BRU 


WRITE 


SKD 


ADIN 


BRU 


WRITE 


BRU 


SEARCH 


ETR 


=0777777 


STA 


*AD0UT 


MP0 


AD0UT 


MP0 


ADIN 


LDA 


*ADIN 


MP0 


ADIN 


LDB 


*ADIN 


LRSD 


6 


STA 


«AD0UT 


MP0 


AD0UT 


LRSB 


6 


STB 


*AD0UT 


MP0 


AD0UT 


LDA 


LASTW 


ETR 


=077777 


SKU 


ADIN 


BRU 


WRITE 


SKG 


ADIN 


BRU 


WRITE 


BRU 


SEARCH 


MP0 


ADIN 


LDB 


*ADIN 


LLSD 


12 


LRSD 


6 


STA 

MP0 


•»AD0UT 

AD0UT 


LRSB 


6 


STB 


*AD0UT 


MP0 


AD0UT 
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MP0 


ADIN 




LDA 


♦ADIN 




LRSA 


6 




STA 


♦AD0UT 




MP0 


AD0UT 




MP0 


ADIN 




LDA 


LASTW 




ETR 


=077777 




SKU 


ADIN 




BRU 


WRITE 




SKG 


ADIN 




BRU 


WRITE 




BRU 


SEARCH 


srchU 


MP0 


ADIN 




LDB 


♦ADIN 




LLSD 


18 




LRSD 


6 




STA 


♦AD0UT 




MP0 


AD0UT 




MP0 


ADIN 




C0FI 


(B,A) 




LDB 


♦ADIN 




LRSA 


12 




LLSD 


6 




STA 


♦AD0OT 




MP0 


AD0UT 




LRSB 


6 




STB 


♦AD0UT 




MP0 


AD0UT 




MP0 


ADIN 




LDA 


LASTW 




ETR 


=077777 




SKU 


ADIN 




BRU 


WRITE 




SKG 


ADIN 




BRU 


WRITE 




BRU 


SEARCH 


IN 


RES 


06000 


0UT 


RES 


06000 


ADIN 


PZE 




AD0UT 


PZE 




LASTW 


PZE 






END 


READ 
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